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We introduce a diffusion modef for energetically inhomogeneous systems. A random walker moves 
on a spin-S Ising configuration, which generates the energy landscape on the lattice through the 
nearest-neighbors interaction. The underlying energetic environment is also made dynamic by prop- 
erly coupling the walker with the spin lattice. In fact, while the walker hops across nearest-neighbor 
sites, it can flip the pertaining spins, realizing a diffusive dynamics for the Ising system. As a result, 
the walk is biased towards high energy regions, namely the boundaries between clusters. Besides, 
the coupling introduced involves, with respect the ordinary diffusion laws, interesting corrections 
depending on either the temperature and the spin magnitude. In particular, they provide a further 
signature of the phase-transition occurring on the magnetic lattice. 
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PACS numbers: 5.50.+q,05.40.Fb,05.60.-k 



I. INTRODUCTION 



The problem of diffusion on inhomogeneous media is attracting much attention, due to its fundamental importance 
in nearly every field of science and engineering [l|, 0, 0, 0, [H, 0] . In fact, for a proper description of many diffusive 
systems (for example proteins "sliding" on DNA, diffusion of charge carriers in solids, flow through porous media, 
etc), apart from the motion of the particles, also the underlying environment must be included [IjQ. 

There are several ways to introduce disorder or, more generally inhomogeneity. For example, it can be geometric, 
due to an irregular lattice structure, or energetic. In the latter case, the lattice sites (or bonds) are assigned different 
energy states and, consequently, the walker is biased towards sites corresponding to potential wells (or small energy 
I , barriers). Moreover, disorder can be deterministic or random and it can be dynamic (the environment is renewed at 
each jump of the walker) [l|, or static (the environment is frozen in a particular configuration) pHIT^, Particles 
diffusing on such structures can also be endowed with memory effects [Ultllj], or be influenced by the distribution of 
other diffusing particles on the same structure [l5j . 

In our model, an evolving, inhomogeneous energy landscape is introduced, by coupling the random walk with a 
spin-S Ising system. More precisely, we assume a spin-S arbitrary lattice and we let a random walker moving on it. The 
relevant energy landscape is then obtained by relating each lattice site with the pertaining nearest-neighbor interaction 
according to the Ising Hamiltonian. In other words, the walker moves on a lattice where each site is occupied by a spin 
Oi G [— S, S] and which generates the energy environment through the Ising interaction. Now, if we make the Ising 
, ferromagnet be in contact with a heat-bath, by varying the temperature parameter the spin configuration evolves 
£N) ' and then, also the energy landscape is modified. In particular, the temperature acts as a dispersion parameter [Tg], 
being able to control the roughness of the energetic environment. In fact, when the temperature is sufficiently low, 
the lattice is ferromagnetic and the energy landscape is flat, viceversa when T — » oo the energy landscape is rugged. 

However, differently from the dynamic, inhomogeneous systems introduced in previous works P, [Io| . where the 
energy landscape was updated from external forces, here we assume that the random walker, while hopping across 
the sites of the underlying lattice, flips the relevant spins. Hence, during the diffusion of the walker on the lattice the 
• i-h , magnetization and the energy properly vary. In fact, as we will see later, by defining a suitable spin-flip probability, the 
random walker is able to provide a diffusive thermal dynamics (l7l [l8|. In particular, as a result of our assumptions, 
the walker is now biased towards such sites that, by flipping the relevant spin, an energy gain can be achieved. 

Therefore, the problem of the RW on an inhomogeneous energy landscape is non trivially extended to the problem 
of their interaction: the RW affects, and is biased by, the energy landscape. In other words, there are two interplaying 
stochastic processes: the motion of the walker and the evolution of the spin configuration. 

Our work will be mainly numerical and the algorithm implemented is very general, being easily appliable to arbitrary 
lattices, made up of spins which can assume an arbitrary, finite number of states. 

The aim of this work is then to characterize the random walk introduced, especially highlighting how its interaction 
with the magnetic lattice affects its diffusion. In particular, it would turn out to be interesting to relate the behavior 
of the walker with the evolution of the energy landscape, namely with the evolution of the magnetic lattice. Hence, 
we analyze our biased random walker (BRW) at different temperatures and then we compare results with those, 
already known, relevant to the ordi nary , unbiased random walker (URW). Interestingly, as we will show, though their 
asymptotic behaviors agree (l9l [20l . |2l| . temperature dependent corrections have to be introduced. In particular, the 
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functional laws describing the behavior of our BRW are URW-consistent, while the pertinent multiplicative factors 
peak at T c . Therefore, effects due to the coupling between the walker and the magnetic system are strongest as 
the latter undergoes its phase transition. In other words, the diffusion of the BRW provides signatures of the phase 
transition occurring on the magnetic lattice. Besides, in order to understand to what extent the spin magnitude 
influences the walker diffusion, we take into account both spin-1/2 and spin-1 Ising systems. 

Finally, notice that all the measures that are being explained are performed after the magnetic system has reached 
a steady state. 

The layout of the paper is as follows. In Sec. [II] we explain how the energy landscape is generated and how our 
BRW can update it; we also underlines the differences with respect to the URW. In Sec. lIIII we show how, under some 
conditions, such differences can vanish and then the URW is recovered. In Sections HVtiVIII we describe the numerical 
simulations performed, useful to characterize the walker behavior. We especially analyze in details the covering time, 
the number of returns to the origin and of distinct sites visited since they better emphasize the relationship between 
the walker and the magnetic lattice. Finally, Sec. I Villi contains a summary and a discussion of results. 



In this work we deal with a RW moving on, and interacting with the energy landscape generated by the following 
Hamiltonian applied to the magnetic configuration of a spin-S Ising system: 



The spin variable a may take the (2S+1) values —S, —S+ 1, . . ., S— 1, S and Aij is the adjacency matrix associated to 
the arbitrary network where spins are placed on. Hence, the first sum only involves nearest neighbor pairs, according 
to the chemical distance. 

Though our analysis has been performed on a toroidal squared lattice with J = 1, h = in order to focus the 
attention on the very dynamical effects, in the remaining of this section we make assumptions on neither the structure 
of the lattice nor on the spin magnitude (though finite). 

Hitherto we have just explained how the energy landscape is generated starting from a discrete spin configuration, 
while now we will describe how the coupling between the magnetic lattice and the walker works. 

The random walker is assumed to be able to move on nearest-neighbor sites or stop, and it can also flip the 
spin pertaining to the reached site (notice that, when the spin magnitude is very large, the latter procedure can 
be quite complex due to a (25* + l)-manifold choice). Therefore, our model displays two interplaying stochastic 
processes: the diffusion of the walker on the lattice and the evolution of the spin configuration. Such processes can be 
considered consequentially (firstly decide the site to move towards and then select the relevant spin state or viceversa) 
or contemporary (consider all possible combinations spin+site and choose one of them). Indeed, in any case, there 
exist many different ways to rule this system, ranging from completely random to completely deterministic. 

The assumptions for our model have been taken in order to realize a proper diffusive dynamics for the Ising model. 
Such a dynamics was introduced in a previous paper [171 ] where the thermodynamics aspects were investigated. In 
particular, it was found that our diffusive thermal dynamics is actually able to drive the system towards a non 
canonical equilibrium state, which depends on the temperature but not on the particular initial spin configuration. 
As far the critical behavior, it preserves the universality class, though the critical temperature is increased: 



Because of this sort of right-shift with respect the canonical dynamics, it is worth underlining that, in the remaining of 
the paper, when we refer to the critical range or temperature, it is always meant according to the diffusive dynamics. 

Now, let us see in detail how the probability running our RW is defined. First of all, it contemporary takes into 
account the motion of the walker and the spin-flip procedure, besides, it is local since it only depends on the magnetic 
configuration of RWs nearest-neighbor sites. More precisely: 



II. DIFFUSIVE DYNAMICS 




(1) 



T S=l/2 = 2.602(1) 
T c s=1 = 1.955(2). 



(2) 
(3) 



'PT(s,j\s ,i) 



PT(s,j)(Ajj +5jj) 



(4) 



represents the probability that the walker, being on site i with coordination number Zi, jumps on a n.n. site j and 
realizes the magnetic configuration s. The spin configuration before the jump is denoted as sq, while {s'} is the set 
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of the new possible configurations. Furthermore, 

PtU^= 1 + Jae As)] (5) 
is derived from the usual Glauber probability (see for more details). Also, 

A25.,(S) = (a} (6) 

is the energy variation consequent to the process, where crj and erj represent the spin-state on site j before and after 
the flip procedure, respectively. 

You can notice that, at each step, the walker can choose among Zi + 1 sites to move towards (or stay on) and, 
contemporary, it can also choose if flip the relevant spin, being biased in order to achieve an energy gain. 

Hence, all in all, there are {zi + 1) x (2S + 1) options including that the magnetic configuration of the system, as 
well as the position of the walker, will possibly remain unchanged. 

Of course, in a (i-dimensional hypercubic lattice, the number of nearest-neighbors does not depend on the particular 
site and Zi = 2d, \fi. 

Notice that, the hopping rate between two sites is, in gen eral, different going forward and backwards, so that the 
random walk is asymmetric. Consequently, as stressed in [l7|, this kind of dynamics does violate the detailed balance 
condition and the equilibrium states achieved are non-canonical. 

In traditional models of diffusion on energetic landscapes, the jump rate is typically controlled by the local energy 
at the start point or by the energy-barrier height between start and end points [12j . Though previous equations 
imply that spin-flips occur on the site where the walker is moving towards, we can as well think our model in terms 
of energetic barriers. In fact, energy-barriers are lower for nearest-neighbor sites which let, by means of spin-flip, 
a higher energy gain. Furthermore, such behavior of the walker is consistent with the physical systems which have 
inspired the model [l7lllg||. 

It is now worth comparing our RW with the traditional unbiased random walker, usually defined according to the 
probability: V(i,j) = — since, at every step the walker must move and the hopping probabilities are isotropic and 
do not depend on time. In this work, in order to establish a stronger analogy with our BRW, we will endow the 
unbiased random walker with a waiting probability so that, from site i, the possible, equivalent, choices are + 1. In 
other words, we allow repetitions within the succession defining the trajectory of the walker and 

Pi'-./; X l ; (?) 

Zi + 1 

where j is a nearest neighbor of i's or, possibly, the site i itself. In the following, we will refer to this able-to-stop 
unbiased random walker as SURW. 

It is known [22j that, for the URW, the possibility of staying on the same site is crucial in the short time regime, 
while in the long time behavior it has no important consequences. As we will see later, an analogous long-time effect 
is also experienced by our BRW. 



III. BRW RECOVERS URW 



As mentioned in Section^ the temperature parameter can tune the roughness of the energy landscape. In particular, 
when T is sufficiently low, the magnetic lattice is homogeneous, the energy landscape is flat and we expect to recover 
the URW case. On the other hand, when T —> oo, a completely disordered lattice and, consequently, a rugged energy 
environment, is achieved. Nevertheless, since the energy variations consequent to whatever possible spin-flip would 
be very small compared with /3, we again expect to recover the URW case. Then, in this section, we want to prove 

N 

that Eq. (|U) recovers Eq. Q, under the conditions (3 — > oo and | cr^ | = NS, or j3 — > 0. Firstly, let us consider the 

i=0 

former case with Oi = S,\/i. Suppose that the walker jumps from site i to j, with coordination numbers Zi and Zj 
respectively and that, consequently, <jj is flipped in a 1 . Then, Eq. (U) can be rewritten as 

P(i,j,<Tj - a') = — , (8) 

EE 

fc=0 ;=o 



FIG. 1: (Color online) Magnetic (left panel) and Visit (right panel) lattice for a 100 x 100 spin-1 Ising system subject to 
the diffusive dynamics described in Sec. [TTJ The sample, initially paramagnetic, was suddenly cooled at low temperatures 
(T < T c ). Note that null spins (white) arrange themselves just on those sites such that their n.n. provide, all in all, a neutral 
magnetization; these sites usually belong to cluster boundaries. In the left panel the warmest colors are reserved to the most 
visited sites which are the ones corresponding to boundaries between clusters; analogous results hold for spin- 1/2. 



where E = e J @/ s and z's nearest-neighbors have been numbered from (the walker remains on i) to 2j. We also 
dropped the factor (Aij + Sij), because we assume j to be linked to i, or, possibly, i — j. Now, since (3 — > oo, then 
E — > oo and we can write: 

V(i ia--a')-i^ + ° iE ~' ) [i<T ' = S (9) 



where 



( = min (z k ) 

k—0,...,Zi 



Conversely, when (3 — > 0, then E — > 1, and, recalling that S is finite, you can easily find that: 



where 



max (z k ). 

k=0,...,Zi 



The previous equation depends neither on a', nor on the magnetic configuration of the lattice and hence, all in all, 
the probability of jumping from a site to another recovers Eq. |7|). 

In the following sections we analyze the behavior of the walker introduced, focusing the attention on those aspects 
which are mostly affected by its interaction with the magnetic lattice. Results will be further stressed by comparison 
with their SURW counterparts. 



IV. VISIT LATTICE 



In this work a task of ours is to relate the motion of the walker with the magnetic configuration of the Ising lattice 
representing the energy landscape. To this aim, we introduce the visit lattice, meant as the L x L array whose 
elements are incremented by a unit each time the walker passes through the pertaining site. In Fig. [1] such a lattice 
is compared with the magnetic one. From Eq. we expect the walker to be attracted towards high energy regions 
which, in our model, corresponds to borders between clusters. Of course, this attraction would affect the distribution 
of visit numbers on the lattice provided that the parameter (3 is not so small to make any spin-flip equally probable 
(see Eq. ([TO])). Actually, in Fig. [TJ the attraction felt by the walker is strong enough to generate detectable effects 
on the visit lattice: as expected, the most "popular sites" are just those belonging to the perimeter of cluster. As a 
consequence, the visit lattice mirrors the magnetic lattice: looking at the former one can derive the spin configuration 
and vice versa. 
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FIG. 2: (ejoc) — (e) for a 240 x 240 spin-1/2 and spin-1 Ising systems initialized ferromagnetic and then heated. The two plots 
can be easily distinguished, as they peak at the relevant critical temperature (dashed line) and the former case shows a higher 
value of T c . 



In the previous section we showed that, according to Eq. [4] BRW does not move freely, but it can be forced to 
stay nearby high energy regions. Hence, we expect the local energy e; oc to be larger than the energy of the whole 
system e. Now, we wonder if the difference between such quantities is somehow temperature dependent. Therefore, 
we consider the quantity e = (e; c) — (e), where, we recall, ti oc represents the energy relevant to the site i occupied 
by the walker, namely ei oc = en a i- As shown in Fig. [2J as long as the temperature is sufficiently low, the lattice 



appears homogeneous and e is null. However, heating the sample, some domains develop and, since the walker verges 
on their borders, (e; oc ) can increase more than (e) so that e rises. While approaching the critical temperature, more 
and more clusters arise and the walker is more and more likely to be found on their boundaries, which explains the 
maximum in T c . Conversely, at high temperature, when the paramagnetic phase has been reached, (e) gets to (e; oc ). 

Note also that, in Fig. [21 the peak relevant to the spin-1 case is sharper, which means a stronger interaction between 
the walker and the magnetic lattice. 

We also measure e for an unbiased random walker allowed to rest and moving on an Ising lattice subject to a 
non-diffusive dynamics. Of course, in this situation, the walker is completely useless for the evolution of the system, 
nevertheless its behavior underlines that results obtained for the BRW are really due to its interaction with the 
magnetic system. In fact, we find that, for the SURW, e remains close to zero without displaying any significant 
dependence on the temperature. 

Therefore e provides a signature that, as far our BRW, at the critical temperature interesting phenomena occur, 
not only in thermodynamics terms. 



The previous two sections pointed out that the hopping-ffipping probability, defined in Section [II actually biases 
the walker towards high energy regions and that, in the critical range, the coupling between the walker and the 
magnetic lattice is even more important. Hence, the phase transition also emerges from the behavior of the walker; 
this interesting feature will be especially taken into account in the following analysis. In particular, we now consider 
the covering time, namely the time (in unit step) taken by the walker to visit all N sites making up the lattice. We 
recall that the lattice is squared and endowed with periodic boundary conditions, so that the walker can actually 
cover an infinite distance on it. 



V. 



LOCAL ENERGY 




VI. COVERING TIME 
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FIG. 3: (Color online) Covering time versus temperature and size of the lattice for the biased random walker moving on a 
spin- 1/2 (left panel) and spin-1 (right panel) Ising system. In both cases the outline depends strongly on the size and the 
temperature affects the covering time just during the critical regime. 
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FIG. 4: Multiplicative constant a versus reduced temperature for the BRW applied to the spin-1/2 (•) and spin-1 (o) Ising 
systems. Both functions peaks at zero, while (xurw (dashed line) is temperature independent. Note that (within the error) 

, 3=1/2 , S=l 
OtSURW S a BRW S a BRW- 



As depicted in Fig. [3l for both spin-1/2 and spin-1 systems, the covering time TJv measured for the BRW increases 
with the size of the lattice and a temperature dependence is also noticeable. In particular, there is an increase in the 
covering time at about the critical temperature, which has been previously measured [ItJ revealing to be fairly larger 
than the canonical one (Eqs. [2] and [3]) 

As far the dependence on the total number of sites N, it is consistent with the logarithmic law: 



T N = aNQogN) 1 



(11) 



found by analytical [23|, |24J as well as numerical [25| methods applied to an unbiased random walk. Then, the very 
effects due to the bias have to be tracked down in the multiplicative factor a. In fact, by fitting, according to Eq. (|ll|l . 

the data relevant to both spin systems, we evidenced that a^j^ 2 and ot^j^ w depend on T and are larger than ajjnw 
found in [23l [25j | . Notice that either the possibility of maintaining the same position and the interaction with the 
magnetic lattice concur in lengthening the covering time, but the role played by the latter is non trivial. In particular, 



FIG. 5: (Color online) Average number of returns to the origin (left panel) and of distinct sites visited (right panel) versus 
temperature and number of steps made by the BRW on a 240 x 240 spin-1 Ising lattice. Analogous results were found for the 
spin-1/2 system, though the critical phenomena are less emphasized. Note that, as n <C N, the walker can not realize the 
finiteness of the lattice. 

it makes and afj^V exhibit a maximum at about T c , namely, in the critical region, it takes more time the 

BRW to cover the lattice. On the other hand, the SURW displays a covering time still consistent with Eq. (fTTj) but, 
of course, independent on T. Note that, as shown in Fig. dj within the error (< 2%), 

, £1=1/2 / S=l ft n\ 

olsurw < at BRW < a B RW ( 12 ) 

As expected, the BRW, with respect the SURW, is slowed down since it may be "trapped" nearby high energy regions 
constituted by sites where several, energetically favorable, sin-flip are possible. On the other hand, the quantities in 
Eq. I|12p are all comparable at low temperatures. In fact, during the ferromagnetic phase, when T <c T c , the lattice 
appears homogeneous and the bias has no effect; an analogous phenomenon is expected at very high temperatures (in 
section Hill we proved that, indeed, in these cases the BRW recovers the SURW). Besides, since in the spin-1 case the 
walker has to manage a greater number of possibilities, the slowing down effect is even higher and the relevant peak 
in Fig. 2] is more marked. 

Note that, with a non-diffusive dynamics, the time required to scan each lattice site at least once can be much 
smaller. For example, by adopting the type-writer sequence, T/v is reduced by a factor Aasnw (logL) 2 . For this 
reason, our diffusive dynamics may be though as "slow" . 

VII. RETURNS TO THE ORIGIN AND DISTINCT SITES VISITED 

In this Section, we want to deal with other two characteristic quantities concerning the BRW: the number of returns 
to the origin R n and of distinct sites visited D n , after an n-step walk. Actually, we ought to distinguish between two 
regimes: 

1. if n is a step number satisfying n <§C N, since the walker will not have sampled a substantial number of sites of 
the lattice, the lattice will appear to be infinite; 

2. in the long time (N <C n) the walker will appreciate the toroidal effect due to the periodic boundary conditions. 

Therefore, according to the walk-length, our results will be compared with those analytically known and relevant to the 
URW on an infinite or periodic 2-dimensional lattice, respectively. As we will see, in both ranges, the exponents found 
for R n and D n agree with known results, while, as far the multiplicative factors, one has to introduce a dependence 
on T. 

Let us firstly consider the short time case. Known results [l!| [2(| about the average number of returns to the origin 
R n and the number of distinct sites visited D n , after an n-step walk for an URW, state that: 



R n - a R logn, 



(13) 



8 



2.1 - 



2.05 



O-OO' 



1.95 



1.9 



1.85 



O 
O _ r, 



1.8 L 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 

Tr 



FIG. 6: Fit parameter ao versus reduced temperature relevant to the BRW on a 240 x 240 spin-1/2 and spin-1 Ising lattice. 
Both functions exhibit a minimum at about the relevant critical temperature (T c m 2.6 and T c ~ 1.96 respectively) which is 
deeper in the latter case. Conversely, at low temperatures, plots overlap. 



D n ~ ™ (14) 
Logn 

where cir and ajy are constant. By fitting and comparing our outcomes relevant to spin-1/2 and spin-1 systems 
(Fig. O, we find that Eqs. IT51 and [Til formally still hold, but an and au are functions of the temperature. More 
precisely, in T c they show a maximum and a minimum, respectively; such effect is more important when S = 1 
(Fig.® . 

When the walk-length is large enough, for the walker, to experience the toroidal effect, we find that the number of 
returns R n recovers the URW case, being 

Rn ~ (15) 

while the number of distinct sites visited is consistent with 

D n ~ N (1 - e- nA °), (16) 

relevant to the URW on a periodic lattice [20], provided that Ad depends on T (Fig. [5]) However, by further increasing 
n, the bias effect vanishes also for D n which, finally, equals N. 

Therefore, our results are asymptotically independent on the temperature and on the spin magnitude, which means 
that the walker looses memory of its bias. 

It should be underlined that, in both regimes, the extreme points recorded at T c for i?„, as well as D n , are consistent 
with what previously found about 7/v and e. 

Analogous results are expected also for the probability of return to the origin Po, n - In fact, we anticipate that the 
BRW still recovers the unbiased exponent, with a multiplicative factor maximum at T c . 

Notice that recovering the conventional diffusive reg ime from a significantly altered model (at least in a long time 
limit) is consistent with several previous works [l|, |26| . The fact that our model yields diffusive behavior should be 
related to the absence of strong memory effects which, indeed, could determine anomalous effects [l3l flij. 



VIII. CONCLUSIONS 



By the analysis performed so far, we are able to characterize the BRW introduced and also to relate its behavior 
with the evolution of the underlying energetic environment. 
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FIG. 7: (Color online) Average number of returns to the origin (left panel) and of distinct sites visited (right panel) versus 
temperature and number of steps n made by the BRW on a 240 x 240 spin-1 Ising lattice. In this case N <C n and the walker 
can experience the periodic boundary conditions the lattice is endowed with. Note that the temperature dependence is now 
scarcely detectable. Similar results were found for the spin- 1/2 system. 
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FIG. 8: Fit parameter Ad versus reduced temperature pertaining to the BRW on a 240 x 240 spin-1/2 and spin-1 Ising lattice. 
Note that the minimum occurs at about the critical temperature (T c w 2.6 and T c ~ 1.96 respectively) and it is deeper in the 
latter case. Conversely, at low temperatures, plots are overlapped. 



Our measures of local energy show that, according to the algorithm introduced in Sec. HH the walker aims to move 
towards high energy regions, where favorable spin-flips can occur. These regions correspond to boundaries between 
clusters, whose concentration depends on the temperature, being largest around T c . Thus, there exist a sort of 
temperature sensitive traps, the walker is attracted to. This attraction reflects on a non-homogeneous visit numbers 
distribution and also causes a slowing down well evinced by measures of covering time, numbers of returns to the 
origin and of distinct sites visited. In general, these quantities asymptotically agree with results, analytically known, 
relevant to the URW. In particular, when the walk-length is not so large to visit all lattice sites, their functional 
forms are URW-consistent, but temperature dependent multiplicative factors have to be introduced. Their particular 
dependence on T has two main consequences: first, the slowing down is not only due to a non-null waiting probability, 
but it is mainly a consequence of the bias introduced; second, the interaction between the walker and the magnetic 
lattice is stronger in the critical region. More precisely, in that temperature range, the length of borders between 



10 



clusters is large so that there are lots of high-energy sites and, contemporary, the temperature is still not too large 
to have a significant (with respect /?) energy gain by spin-flips. Hence, multiplicative factors just peak at the critical 
temperature. 

What has been said hitherto holds for either spin-1/2 and spin-1 systems. In fact, our analysis have been contem- 
porary performed on both, in order to evidence how the spin magnitude affects the walker diffusion on the lattice. 
Then, our results show that peaks get sharper and higher when the number of spin states is larger. Actually, in the 
latter case, the walker has to manage with more opportunities and, consequently, it is further slowed down. 

In summary, in the limit T — > (T — > oo) the underlying energetic landscape becomes homogeneous and the BRW 
recovers the case of an ordinary unbiased random walker endowed with a non-null waiting probability. On the other 
hand, when T = T c we record the most important effects. In fact, T c is an extremal point for the multiplicative factors 
which correct the ordinary laws. Hence, the bias introduced determines stronger effects as the critical temperature is 
approached, though not affecting the diffusive regime. Therefore the BRW behavior gives a further evidence of the 
phase transition. 

Since such corrections are more important when the spin magnitude is larger, we argue that an investigation in the 
continuum limit for S would turn out to be useful in order to clear the nature of the above mentioned extreme points. 



I. Avramov, A. Milchev and P. Argyrakis, Phys. Rev. E 47, 2303 (1993). 

J. Bricmont and A. Kupiainen, Phys. Rev. Lett. 66, 1689 (1991). 

A. Nitzan and M. A. Ratner, J. Phys. Chem. 98, 1765 (1994). 

S. Bustingorry, E. R. Reyes, and M. O. Caceres, Phys. Rev. E 62, 7664 (2000). 

A. K. Harrison and R. Zwanzig, Phys. Rev. A 32, 1072 (1985). 

S.D. Druger, M.A. Ratner and A. Nitzan, Phys. Rev. B 31, 3939 (1985). 

N. Masuda and N. Konno, Phys. Rev. E 69, (2004) 66113. 

M. Slutsky, M. Kardar and L.A. Mirny, Phys. Rev. E 69, 61903 (2004). 

B. D. Bookout and RE. Parris, Phys. Rev. Lett. 71, 16 (1993). 
LA. Levitsky, Phys. Rev. B 49, 15594 (1994). 

E. Arapaki, P. Argyrakis, I. Avramov and A. Milchev, Phys. Rev. E, 56, 56 (1997). 

W. Haus and K.W. Kehr, Phys. Rep. 150, 263 (1987). 

B. Schulz, S. Trimper and M. Schulz, Eur. Phys. J. B 15, 499 (2000). 

B. Schulz, S. Trimper, Phys. Lett. A 256, 266 (1999). 

M. Schulz, S. Trimper, Phys. Rev. E 62, 221 (2000). 

K. Sapag, V. Pereyra, J. L. Riccardo and G. Zgrablich, Surf. Sci. 295, 433 (1993). 
E. Agliari, R. Burioni, D. Cassi and A. Vezzani, Eur. Phys. J. B, 46, 109 (2005). 

P. Buonsante, R. Burioni, D. Cassi and A. Vezzani, Phys. Rev. E, 66, 36121 (2002), and references therein. 
E.W. Montroll and G.H. Weiss, J. Math. Phys. 6, 167 (1965). 

G. H. Weiss, Aspects and Applications of the Random Walk (elsevier Science, Amsterdam 1994) 

H. Weiss and R.J. Rubin, Adv. Chem. Phys. 52, 363 (1983). 
R. Burioni and D. Cassi, J. Phys. A 38, (2005) R45. 

J. Jonasson and O. Schramm, ArXiv Mathematics e-prints, arXiv:math/0002034 (2000). 

A. Dembo, Y. Peres, J. Rosen and O. Zeitouni, ArXiv Mathematics e-prints, arXiv:math/0107191 (2001). 

A.M. Nemirovsky, H.O. Martin and M.D. Coutinho-Filho, Phys. Rev. A 41, 761 (1990). 

S. Trimper, U.C. Tauber, and G.M. Schiitz, Phys. Rev. E 62, 6071 (2000). 



